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^ I Abstract 

"q ■ In this article the stability of the Standard Model (SM) vacuum in the presence of 

radiative corrections and for a Higgs boson with a mass in the vicinity of 125 GeV 
■^j- \ is discussed. The central piece in this discussion will be the Higgs self-interaction A 

and its evolution with the energy scale of a given physical process. This is described 
^ ■ by the /3-function to which we recently computed analytically the dominant three-loop 

^ contributions [I] J! These are mainly the QCD and top- Yukawa corrections as well as 

Ok the contributions from the Higgs self-interaction itself. We will see that for a Higgs 

boson with a mass of about 125 GeV the question whether the SM vacuum is stable 
and therefore whether the SM could be valid up to Planck scale cannot be answered 
CN " with certainty due to large experimental uncertainties, mainly in the top quark mass. 
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1 The Higgs potential and the stability of the SM 
vacuum 

o 



In the SM the Higgs potential at tree-level appears as part of the Lagrangian for a 
scalar SU(2)-doublet field: 

C* = d»&& l $- (mV^ A (V<I>) 2 ), $=^| 1 ^J. (1) 

This doublet aquires a non-zero vacuum expectation value (VEV) under spontaneous 
symmetry breaking (SSB) and we get the Higgs field, three Goldstone bosong^l and the 



1 email: max.zoller@kit.edu 

2 In pQ we also give the dominant contributions to the /3-functions for the top- Yukawa coupling, the 
strong coupling and the anomalous dimensions of the scalar, gluon and quark fields in the unbroken 
phase of the Standard Model at three-loop level. 

3 These Goldstone bosons can be absorbed by the massive W and Z bosons. 
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masses of the SM particles: 
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At tree level the mass of the Higgs boson is then given by 

M 2 H = -2m 2 = 2Xv 2 . (3) 
For Mjj = 125 GeV the Higgs potential is shown in Fig. [TJ If we now include radia- 
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Figure 1: The SM Higgs potential at tree-level for Mh = 125 GeV 

tive corrections we have to consider the effective potential V r e //(A(A), ^(A), $(A)) as 
introduced in [2]. All couplings and fields undergow an evolution up to some scale A 
where the theory seizes to be validfl Here we want to investigate whether a scenario 
in which the SM provides a good description of nature up to the Planck scale, i.e. 
A ~ 10 18 GeV, is possible!! The general shape of this effective potential for the SM is 



M H >m min 



M H <m mill 



Figure 2: The effective Higgs potential 
shown in Fig. [2] for the cases of a Higgs mass larger and smaller than a critical value 



4 With t = log (j^j we have the held $(A) = $ c/ -exp (f dt'^(X(t'), gi(t'))dt'j where $ c ; is the 

classical held in the absence of radiative corrections. fiQ is the scale at which we start the running of 
the paramters. 

5 An extended model is expected to be needed at this scale due to gravity. 
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m m i n , the minimal stability bounc{§ (see also [3]). If the Higgs mass is chosen below 
this critical value a second minimum develops which is lower than the SM one. This 
implies that the SM vacuum is no longer stable, i.e. it can tunnel into this energeti- 
cally favoured state. Since this is in contradiction to our observation^ we are led to the 
conclusion that our theory is incomplete and that new physics has to enter between 
Fermi and Planck scale. It has been demonstrated in [I] that for $ ~ A ^> /i a good 
approximation for the effective potential is 

y e// [$]^A(A)$ 4 + 0(A 2 (A),^(A)), (4) 

which means that the stability of the SM vacuum is approximately equivalent to the 
question whether A stays positive up to the scale A (see also [S1E])- Fig- El shows the 
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Figure 3: The behaviour of A(/i) for different Higgs mass values 



evolution of A for different Higgs mass values. For a large enough Higgs mass A diverges 
quickljj§ and for a small Higgs mass A would become negative at a relatively low scale 
before eventually reaching a Landau polej^l The interesting region is around 125 GeV 
where A is very close to zero at the Planck scale! As the minimal stability bound has 
been estimated to be about 129 ± 3 GeV [TO], which is very close to the mass of the 



6 There is also an upper bound m max on the Higgs mass from the requirement that no Landau 
pole appears at energies /i < A. 

7 It is however not possible to exclude that we live in metastable universe, i.e. the lifetime of the 
(local) SM minimum could be longer than the age of the universe. 

8 We find a Landau pole below the Planck scale for Mh > m max ss 175 GeV (see also [S1[7HS])- 

9 The reason for the fact that A becomes positive again and diverges is the evolution of the gauge 
couplings and y t . At large scales (above 10 16 GeV) the electroweak couplings, especially the U(l) 
coupling <7i, start to be the dominant contributions. In contrast y t , which is responsible for the 
decrease of A at lower scales, becomes small. This is also the reason why the effective potential goes 
up again at large field strengths $ beyond the second minimum. 
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boson recently descovered at the LHC, the question of vacuum stability in the SM 
becomes one for precision physics. This serves as a strong motivation for calculating 
the three-loop /3-function of the Higgs self-interaction which describes the evolution of 
this crucial coupling as well as the /3-functions of the relevant SM parameters on which 
j3\ depends. 



2 The three-loop /3-function for the Higgs self- 
interaction 

The /3-function is defined as the derivative of the coupling with respect to the renor- 
malization scale: 

P\(\y t ,9i,---) = ^j-^Kv)- ( 5 ) 

This object has been known at the two-loop level for a while including the dependence 
on all the gauge couplings gi,g2 and g s , the quartic Higgs self-interaction A and the 
Yukawa couplings which give mass to the fermions [TTHT3] . Recently, the /3-functions 
for the gauge couplings have been computed at three- loop level [H]. At two-loop level 
they had been been derived in [TBTfTS] . As there are many Feynman diagrams at the 
three-loop level and the treatment of 75 matrices in dimensional regularization poses 
a serious problem at this order for the Higgs and Yukawa sector we try to find the 
dominant contributions to the evolution of the Higgs self-interaction at the scales that 
we are intersted in. At the scale of the top mass \x = M t = 172.9 GeV we find the strong 
coupling g a pa 1.17 and the top- Yukawa coupling y t pa 0.93 to be much larger than the 
electroweak couplings g^ ~ 0.65 and g\ pa 0.36. The Higgs self-interaction for a Higgs 
boson around 125 GeV is \{M H ) pa 0.13. From this we conclude that a simplified 
model containing only g s , y t and A will give the numerically largest terms for (3\ at 
three-loop level. When combinen with the full SM result at one and two-loop level 
this leads to the evolution of A to the highest precision so far. The Lagrangian of our 
model consists of three pieces: 

£ = £ QCD + £ ti + £ yt . (6) 

with the standard QCD Lagrangian, the Higgs part as defined in (TjQ) and the top- 
Yukawa sector 

4* = -y t {t R (* 2j -$o ■ [I ) + (t, l) L • ( *j, ) * a j (7) 

The indices L and R indicate the left- and right-handed part of the Dirac fermion fields 
as obtained by the projectors 

^ = ^(1-75) P fl = ^(l + 75 ). (8) 

2 2 

10 These couplings enter at every loop order with characteristic factors f 2 - ~ 0.11, I 2 - ~ 0.03, 
2 2 ^ 

w 0.01 and j 1 - ~ 0.07. The Higgs self-coupling is already present at tree level and enters at 

each order linearly as s» 0.01. The second-largest Yukawa coupling to be considered would be 

Ub = ~ 0.02 which is negligible in comparison. 
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Hence 75 enters into our calculation and has to be treated carefully as described below. 



2.1 Calculation 



In order to obtain the /3-functions for a coupling, here A or y t , we have to compute the 
renormalization constant for a vertex involving this coupling and the field renormal- 
ization constants for the external legs of this vertex. The latter is done by computing 
loop corrections to the respective propagators. The renormalization constant for y t for 
example can be computed as^ 



^(Yukawa vertex) 



~(top field) ^(top field) (Higgs field) 
L H 



(9) 



The /3-function f3 yt can then be computed from the requirement that the bare Yukawa 
coupling 

y b r = z yt [ gi ^)M^)M]y^) (10) 

is independent of the renormalization scale \i and hence its /u 2 -derivative must vanish. 
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Figure 4: Sample Feynman diagrams for the calculation of renormalization constants 
for the quartic Higgs vertex, the Higgs self-energy and a Yukawa vertex at three-loop 
order. 



For the purpose of calculating renormalization constants and /3-functions we only need 
the UV-divergent part of all these diagrams. Nevertheless, there are two issues which 

11 Note that in general left and right-handed fields have to be renormalized separately. 
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have to be considered carefully in a computation of this type. The first problem is the 
complicated nature of our integrals when we have arbitrary momenta flowing into the 
external legs of our Feynman diagrams. A nice feature of the MS-scheme however is that 
renormalization constants do not depend on external moment a[^l. so we can set those to 
zero. Unfortunately, this introduces artificial IR-divergences into our diagrams which 
cannot be distinguished from the UV ones in dimensional regularization. In many cases 
this problem can be avoided by setting all external momenta to zero except for one 
which enters at one leg and exits at another. If we have no masses, as we do here, the 
propagator-like integrals resulting from this method can be computed with the FORM 
package MINCER [T3] up to three-loop order. For the A-vertex and its radiative 
corrections this fails, however, as there are still IR-divergent diagrams. A method to 
compute only the UV-divergences without having to worry about the IR ones has been 
described in pOll^Tj . The trick is to introduce the same auxiliary mass into every 
propagator denominator, Taylor expand in all the external momenta and introduce all 
possible counterterms for the auxiliary mass in order to cancel sub divergences arising 
from this new mass. The resulting massive tadpole integrals can be computed with 
the FORM package MATAD [22] • This method will yield the correct UV-pole part of 
the calculated Feynman diagrams (but not the correct finite part) which is enough for 
the computation of /3-functions. Where possible both methods have been used for the 
calculation in order to have an independent check. 

The second problem is the treatment of 75 matrices, which arise here from the projectors 
in the Yukawa sector, in d space-time dimensions. This matrix is only well-defined in 
four dimensions: 

75 = ijVlV = ^'WtVtV where £0123 = 1 = -£° 123 - (11) 

A naive treatment of 75, i.e. using the relations {75,7^} = and 75 = 1 to eliminate 
as many 75 as possible and then discarding all terms which still have one 75 in them, 
can only be applied to external fermion lines and closed fermion loops with less than 
four Lorentz indices or momenta flowing in or out of the fermion loop in question!^! An 
example for a problematic diagram is shown in Fig. H] (3c). The Lorentz indices of the 
two gluons connected to the closed fermion loop and the two loop momenta pi and P2 
offer the possibility of a non-trivial 75 contribution from such a diagram. And indeed, 
if we apply the treatment suggested by 't Hooft and Veltman in [23], i-e. using the 
definition with e M ^ pCT from (ITTj) and contracting the e-tensor from the closed fermion 
loop with the one from the external fermion line or a projector oc 75 acting on this 
external line to make the integral scalar, we find a sizable contribution. But even this 
treatment is not exact, in fact it is only correct up to an error of 0(e). Fortunately, 
there are only first order poles in e from these integrals which makes the UV-divergent 

12 To be precise, UV-divergent terms can only depend polynomially on external momenta (and 
masses). 

13 Four different Lorentz indices are needed to support a non-vanishing £ M „ p(T after the trace over 
the fermion line has been performed. These can be the Lorentz indices of gauge boson vertices or 
they can be contracted with internal momenta from other loops. External momenta can be set to zero 
without changing the UV-divergent part and can hence be ignored in this consideration. 
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part and therefore our renormalization constants correct. Once again, the finite part 
calculated for these diagrams is unreliable but lukily not needed here. 



2.2 Result 

The result for the /3-function of the Higgs self-interaction 

a 00 1 

M ! ^A(/i) = /3 x (g.,y t ,\) = ^_j-^W(y.,y t ,A) (12) 

in our simplified version of the Standard Model is given by 
/?«=12A 2 + 6y t 2 A-3 ?/t 4 , 

^ = - 156 A 3 - 72 y t 2 A 2 - 2 y ( 4 A + 15 y\ + 40 g^X - 16 
/?f =A 4 (3588 + 2OI6C3) + 873 y*\ 3 + y?\ 2 ( ™ + 756Cs) 

+ rfA(^-19 86 )-rf(if + 366) 

+ s s y A 2 (-1224 + 1152Cs) + tfyfX (895 - 1296C 3 ) 

/ 1 820 

+ 9 y t (-38 + 240C 3 ) + g A a y 2 t \ — - 32n f - 48C; 



V 3 



+ ^ t 4 (-^ + 20n / + 32C 3 ) 



in the MS-scheme. f3^\ Px and (3^ are the one, two and three-loop results re- 
spectively. To get an idea of the size of the individual terms and the overall result 
we evaluate f3 x at the scale \x = M z (with an assumed Higgs mass of 125 GeV and 
the number of fermion flavours n f = 6) which yields a value of j3 x ~ (—0.01) at the 
one-loop level. The two and three-loop contributions change this result by ~ 1% and 
~ (—0.04)% respectively which is quite small e.g. in comparison with the /3-function 
for the top- Yukawa coupling where we get a value of f3 yt ~ (—0.023) and corrections of 
~ 16.6% and ~ 0.7% at two and three-loop level. The full result for j3 yt can be found 
in [TJ . If we have a look at the the numerically largest individual terms in (12. 2|) we find 
a curious behaviour: 

B {3) 

1 r (/i = M z ) = (+7.9 -4.8 -3.1 -2.5 +2.6) x 10~ 5 . (13) 



(16+>) ; 

9 2 a yf yf g 2 yf* giyf giy 2 * 

There is a strong cancellation between these terms making the overall effect almost 
two orders of magnitude smaller than the largest individual contributions This 
significantly improves the convergence of the perturbation series for f3\ and makes the 
remaining theoretical uncertainty small. 



Note that for a significantly different Higgs mass this would not happen to such extent. 
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3 The evolution of the Higgs self-coupling 



Now we want to investigate the effect of the new three-loop result /3i on the running 
of A and therefore its effect on the stability of the electroweak vacuum in the SM. For 
this we also include the electroweak contributions up to the two-loop level. For recent 
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Figure 5: Evolution of A with the scale /j: 2 loop (dashed, blue) and 3 loop (continuous, 
red) results; Uncertainties with respect to the two-loop result: ±la aa , ±l0M t (dotted) 
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detailed discussions of the SM vacuum stability and its dependence on key parameters 
like the Higgs and the top mass see for example [TJEIE1H25]- We start the evolution 
of A at the scale of the top mass and go up to the Planck scale at 10 18 GeV. To find 
the starting values in the MS-scheme we have to match the physical parameters, like 
e.g. pole masses, to their MS counterparts!^! These matching relations depend on the 
exact values of the MS parameter a a {Mz)i the pole mass M t of the top quark and of 
course the mass of the Higgs boson Mh- We use [30J 

a s (M z ) = 0.1184 ±0.0007, M t = 172.9 ± 0.6 ± 0.9 GeV. (14) 

Fig. shows the evolution of A for the two cases Mh = 124 GeV and Mh = 126 GeV.. 
The dependence of the A- running on the parameters a a (Mz) and M t can be esti- 
mated from the shifted curves where these parameters are changed by ±cr as given 
in eq. f ll4p . Note that there is a considerable difference between Mh = 124 GeV and 
Mh = 126 GeV which means that the evolution of A is very sensitive to the value of 
the Higgs mass. Given a fixed value for Mh the largest uncertainty lies in the exact 
value of the top mass. The second largest uncertainty comes from ot a . The total effect 
due to the three-loop part of the /3-functions, an extension of the vacuum stability to 
larger scales, is a little smaller than the a s uncertainty. 

In conclusion it can be said that for a Higgs mass in the vicinity of 125 GeV the 
question of vacuum stability cannot be resolved with certainty. Although it looks as 
if A becomes indeed negative at high scales the experimental uncertainties, mainly on 
the Higgs mass and the top mass, are too large at the moment. However, this is a very 
good motivation to attempt higher precision experiments and calculations in the near 
future. 
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